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Abstract 

We report the value of Og obtained from QCD analysis of existing 
data on deep-inelastic scattering of charged leptons off proton and deu¬ 
terium and estimate its theoretical uncertainties with particular atten¬ 
tion paid to impact of the high-twist contribution to the deep-inelastic- 
scattering structure functions. Taking into account the major uncer¬ 
tainties the value = 0.1143±0.0014(exp.) ±0.0013(theor.) 

is obtained. An extrapolation of the LO-NLO-NNLO results to the 
higher orders makes it possible to estimate ~ 0.113. 
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1 Introduction 


The value of strong coupling constant Ofs can be extracted from the existing 
data on deep-inelastic-scattering (DIS) with experimental uncertainty at the 
one-percent level. With such experimental accuracy achieved the theoretical 
uncertainties become dominant. One of the major theoretical uncertainty in 
as derived from the next-to-leading-order (NLO) QCD analysis of DIS data 
comes from the higher-order (HO) corrections. Therefore, consideration of 
the next-to-next-to-leading (NNLO) corrections is of great importance for 
suppression of the errors in ag. Indeed, account of the NNLO terms in the 
QCD evolution equations for the DIS structure functions, which have been 
calculated recently 0121. allows to reduce the HO uncertainty in ag by a 
factor of 2 -h 3, as it was estimated in Ref. [3]. Another important source 
of the uncertainty in Og is the high-twist (HT) contribution to the structure 
functions. This contribution inevitably appear in the operator product ex¬ 
pansion j3], but unfortunately cannot be reliably estimated hitherto. In the 
limited range of the momentum transfer Q the HT terms, which fall with Q 
by the power law, can simulate the logarithmic-type behavior of the leading 
twist (LT) term. For this reason the magnitude of HT terms is correlated 
with the value of ag, which dehnes the slope of the LT term, and the uncer¬ 
tainty in HT terms propagates into the uncertainty of Og. 

We study impact of these sources of uncertainties on the value of ttg ob¬ 
tained from a NNLO QCD analysis of existing data on the charged-leptons 
DIS off the proton and deuterium targets j3]. The data outside the kine- 
matical region of = 2.5 -i- 300 GeV^ and x = 10“^ -j- 0.75 are left out 
to suppress potentially dangerous theoretical uncertainties keeping sufficient 
sensitivity of data to the value of Og. The data are analyzed using a numer¬ 
ical integration of the QCD evolution equations in the x-space (details are 
described elsewhere jS]). The value of a* obtained in the recent update of 
this £t is [7] 

af^LO(Mz) = 0.1143 ± 0.0014(exp.), (1) 

to be compared with 

af^^°(Mz) = 0.1166 ± 0.0009(exp.) (2) 

obtained by Santiago and Yndurain from the analysis of the proton subset of 
data jH]. The achieved experimental error in a* is quite small and therefore, 
complimentary analysis of theoretical uncertainties is necessary in order to 
clarify the conhdent range of Og variation. 
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Figure 1: The central values of as{Mz) and the experimental errors obtained 
in different orders of pQCD (full symbols). The open symbol shows the result 
of extrapolation to the N^LO. 

2 Uncertainty due to the higher-order cor¬ 
rections 

Generally, estimate of the HO uncertainty is impossible before the HO cor¬ 
rections are known. However, the approach based on variation of the QCD 
renormalization scale (RS) is commonly used for estimation of the HO un¬ 
certainties. This approach is very approximate since it does not account 
for a possible a:-dependence of the corresponding scale-factor and since the 
variation range of the scale-factor is optional. Nevertheless, following this 
approach we estimate the HO uncertainty as the change of Os caused by 
variation of the QCD RS from Q to 2Q. The estimates obtained are given in 
Table |21 Note the reduction of Og with the pQCD order approximately coin¬ 
cide with the corresponding HO uncertainties^. Extrapolating this regularity 

^The decrease in as because of going from the NLO to the NNLO disagrees with the 
results of NNLO analysis by Santiago and Yndurain, but agrees with the estimates of 
Ref.|3] and results of the NNLO fit of Ref.jni. 
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to the N^LO we have an estimate ~ 0.113. 

Table 1: The values of a^iMz) obtained in different orders of pQCD analysis 
and their RS errors. 


LO 

0.1301 ± 0.0026(exp.) ± 0.0149(RS) 

NLO 

0.1171 ± 0.0015(exp.) ± 0.0033(RS) 

NNLO 

0.1143 ± 0.0014(exp.) ± 0.0009(RS) 


3 Uncertainty due to the high-twist contri¬ 
bution 

The HT terms can be taken into account using different approaches. Firstly, 
one can parameterize the HT terms by a flexible model independent function. 
In this approach possible uncertainties in variation of the HT terms allowed 
by data merge into the total experimental error in Og. Secondly, one can cut 
data with low momentum transfer Q and/or hadronic mass W. Evidently in 
the second approach the HT error vanishes, however the statistical error can 
rise because the QCD evolution introduces the largest effect at small Q and 
therefore this region is most sensitive to the value of ttg- Thirdly, one can 
take into account the HT terms using available theoretical models. In this 
approach the HT uncertainty is subject of belief in reliability of the model. 

In the analysis of Ref. [ 7 ] we parameterize the twist-4 terms in additive 
form, F 2 ,l = + if 2 ,L(a^)/Q^, where 

rected for the target mass effects ug. The functions if 2 ,L are dehned in a 
piece-linear form with htted coefficients that provides the model-independent 
account of the HT terms. One does need to include the HT terms in the anal¬ 
ysis since they improve quality of the £t {x^/NDP = 2521/2274 with and 
X^/NDP = 2851/2274 without them). Possible logarithmic Q-dependence 
of H{x) at a; ^ 1 could decrease the htted Og value jTTj, however data used 
in our analysis are not sensitive to such modihcations of H[x) because of the 
cut on X. The value of Og is strongly anti-correlated with H 2 at large x (see 
Figigi. For this reason the error in ag is much larger than it could be in the 
absence of the HT terms. Indeed, in our trial £t with the HT terms set to 
0, the error in ag(Mz) is 0.0003. Comparing this error with the error in Og 
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Figure 2: Shown is the correlation coefficient for the htted values of as{Mz) 
and H 2 at different x. 

obtained in our analysis with the HT terms accounted for, we estimate the 
HT error in ag(Mz) to be 0.0014. This is the major source of the error in 
^NNLO^ more important than the estimated uncertainty due to the HO. 

The results of scan hts with different cuts on Q and with no HT terms 
included are given in Figl^l These results suggest the optimal value of Qmin ~ 
20 GeV^, which provide balance between stability of the central value of 
tts and its error. For Q^in = 21.5 GeV^ the value of as{Mz) = 0.1106 ± 
0.0020 (exp.) was obtained. This value agrees with the results of analysis 
with the HT terms included that supports the both approaches. At that the 
value of tts derived from the £t with the HT terms included is preferable since 
it is more robust due to the wider set of data was used in this case. The cut 
on W leaves out the data at small Q and large x and therefore makes fit less 
sensitive to the HT terms however not so efficiently as the cut on Q since 
the HT terms are important for moderate x as well. With the increase of the 
minimal value of W the error in as gets inappreciably large earlier than its 
central value is sustained. 

In the analysis by Santiago and Yndurain the value of Og was obtained 
from the £t with the HT terms set to 0. The uncertainty in as caused by this 
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Figure 3: The values of as{Mz) obtained from the £t with if 2 ,L = 0 and 
different cuts Qmin (left panel) and the corresponding values of DP 
(right panel). The horizontal lines in the left panel correspond to the central 
value and the error band of as(Mz) in the analysis with the HT parameterized 
in the piece-linear form, the error bars in the right panel are ^J2/NDP. 

constraint was estimated as a difference between the value of as obtained 
from the fit with the HT terms set to 0 and that of the fit with HT terms 
described by the gradient model. In this model the HT contribution to iV 
reads A 2 [dF,^^{x)/dx]/Q^, where is the non-singlet structure function 
and A 2 is a normalization parameter determined from data. 

We hnd A 2 = —0.018 ± 0.006 GeV^ in the trail fit with H 2 {x) described 
by the gradient modeF. We find that in this £t as(Mz) changes by 0.0006 
with respect to that of the £t with HT terms set to 0. This value agrees with 
0.0004 obtained in the analysis by Santiago and Yndurain. However this shift 
in ttg cannot be considered as a reliable estimate of the HT uncertainty in as, 
since the gradient model describes the HT terms very poorly {x‘^/NDP = 
2543/2067). This is not surprising, because the gradient model can only be 

^In the trial fit we leave out the data points with < 3.5 GeV^ in order to better 
match the data set used in Ref.|H]. 
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Figure 4: One-standard-deviation band for the HT term in the proton struc¬ 
ture function F 2 obtained in the gradient model £t (curves) compared to the 
model independent determination (points). 

justihed at x —> 1 and cannot be applied in the whole region of x (see FigE)). 
A correct estimate of the HT uncertainty in ag in the analysis by Santiago 
and Yndurain would be given by a variation of as{Mz) due to going from the 
£t with the HT terms set to 0 to the £t with the HT terms parameterized in 
the model-independent way. In our analysis such variation is —0.0069, much 
larger than estimate of the HT uncertainty based on the gradient model. 

Taking into account this variation, we conclude that in fact Eqn. (CD 
disagrees with Eqn. (ED This disagreement is evidently connected with the 
substantial difference between the gluon distributions obtained in different 
hts. The gluon distribution was estimated in the analysis by Santiago and 
Yndurain as ~ — x)®'^ ^2] with the momentum carried by gluons 

0.752 ± 0.014 at = 12 GeV^. This is much larger than the typical value 
of ~ 0.4 obtained in the recent global hts of PDFs HI eg. Independent es¬ 
timation of this quantity based on the results of CDHS collaboration m is 
0.56 ± 0.02, also much lower than the result by Santiago and Yndurain. We 
comment that due to the value of Og is anti-correlated with the gluon distri¬ 
bution at small x, the increase in Og caused by the absence of the HT terms 
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in that analysis is compensated by its decrease due to enhanced gluons. 


4 Conclusion 

We also estimate theoretical errors connected with account of the heavy 
quarks contribution. These errors are more important than in our previous 
determination of as due to growing impact of the recent low-x HERA data. 
The uncertainty in as{Mz) due to choice between the hxed-flavor-number and 
the variable-flavor-number factorization schemes is estimated as 0.0007. The 
variation of the c-quark mass by 0.25 GeV leads to the variation in as{Mz) of 
0.0006. Other possible sources of theoretical errors, including uncertainties 
in the NNLO evolution kernel and in the strange quarks distribution, are not 
considered because they give much smaller effect. Summarizing all sources 
of theoretical errors we get 

= 0.1143 ± 0.0014(exp.) ± 0.0013(theor.). 

This value is somewhat lower than the world average value of Ref. na but 
agrees with it within errors. 
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